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standard model background expectation and a model independent upper limit on the 
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1 Introduction 

In July 2012 the ATLAS and CMS collaborations at the CERN LHC announced the discovery 
of a particle ICEHSI with properties consistent with the standard model (SM) Higgs boson 0H3- 
Direct measurements of the production and decay rates of the new particle, using SM decay 
channels, have so far played a key role in determining whether or not it is indeed consistent 
with the SM predictions. However, substantially increasing the precision of these measure¬ 
ments will require further data. Searches for Higgs bosons through production mechanisms 
not predicted by the SM, or decay modes involving particles not included in the SM, provide 
a complementary approach and have the advantage of probing specific types of new physics 
models with the existing data. 

This letter presents a search for the pair production of new light bosons (denoted as 'a') de¬ 
caying to pairs of isolated, oppositely charged muons (dimuons). One production mechanism 
for these new bosons is in the decay chain of a Higgs boson h, which can be SM-like or not: 
h —> 2a + X —> + X, where X denotes possible additional particles from cascade decays of 
the Higgs boson. A range of new physics scenarios predict this decay topology, including the 
next-to-minimal supersymmetric standard model (NMSSM) §3 and models with hidden (or 
dark) sectors li^HTlTl. 

The NMSSM is an extension of the minimal supersymmetric standard model (MSSM) Il2l il3fl 
that includes an additional gauge singlet field. It resolves the so-called }i problem Ifl4ll and 
significantly reduces the amount of fine tuning required in the MSSM lH5l . The NMSSM Higgs 
sector consists of three CP-even neutral Higgs bosons hi 2 , 3 , two CP-odd neutral Higgs bosons 
ai ,2 and a pair of charged Higgs bosons h*. The hi and I 12 can decay via h 1 ^ — > 2a\, where 
either the hi or I 12 can be the boson observed at 125 GeV. The ai boson can be light and couple 
weakly to SM particles with a coupling to fermions proportional to the fermion mass. Therefore 
it can have a substantial branching fraction B(a\ — > }i + y.~) if its mass is within the range 
2mp < m a | < 2m t lfl6llT7| (benchmark model 1 in this letter). A search for final states containing 
muon pairs provides sensitivity to models of this form. 

Supersymmetry (SUSY) models with dark sectors (dark SUSY) offer an explanation for the 
excess in the ratio of the positron flux to the combined flux of positrons and electrons observed 
by the satellite experiments lfL8ll20ll in primary cosmic rays as well as predict cold dark matter 
with a scale of 0(1 TeV). A simple realization of these models includes a new U(1)d symmetry 
(the subscript "D" stands for "Dark") which is broken and gives rise to massive dark photons 
(denoted as 70 ). Kinetic mixing of the new U(1)d with the SM hypercharge U(l)y provides a 
small mixing between jo and the SM photon which allows 7 d to decay to SM particles ED . 
Depending on the value e of the kinetic mixing, the 7 d may also be long-lived. The lack of 
an antiproton to proton ratio excess of the magnitude similar to the positron excess in the 
measurements of the cosmic ray spectrum constrains the mass of 70 to be less than twice the 
mass of the proton Ii22l . If the hidden sector directly or indirectly interacts with the Higgs 
field, a number of possible scenarios may be realized. One such scenario, denoted in this letter 
as benchmark model 2, is a model of SUSY where the SM-like Higgs boson can decay via 
h —> 2 nj, where n ] is the lightest neutralino in the visible (as opposed to hidden) part of the 
SUSY spectrum. The ni can decay via ni —» np + 7 d, where np is a dark neutralino that 
escapes detection. Assuming that 70 can only decay to SM particles, the branching fraction 
£?( 7 d — > }i l }i~) can be as large as 45%, depending on the mass of 7 d IflTI . 

Previous searches for pair production of new light bosons decaying into dimuons were per¬ 
formed at the Tevatron [23] and the LHC [24,125]. Searches for associated production of the 
light CP-odd scalar bosons have been performed at e 1 e colliders |26> 27] and the Tevatron H28l . 
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3 Data selection 


Direct a\ production has been studied at the LHC l29l , but this is heavily suppressed by the typ¬ 
ically very weak couplings of the new bosons to SM particles. The constraints on the allowed 
NMSSM parameter space are driven by the measurements of relic density by WMAP H30H and 
more recently by PLANCK (311 , while specifically for the Higgs sector the most relevant mea¬ 
surements come from LEP (32)4371 , LHC measurements of the SM-like Higgs properties, and 
direct searches for h —> aa (25l . In the framework of dark SUSY, experimental searches for dark 
photons have focused on their production at the end of SUSY cascades at the Tevatron (38)4401 
and the LHC (41)1421. Searches at a range of low energy e - e colliders ( KLOE HU, BaBar HU), 
heavy-ion colliders (PHENIX 05)), fixed-target experiments (APEX (461 , A1 at MAMI HZL 
HADES HU), as well as cosmological measurements HU-EU and others (52)4561 provide con¬ 
straints on complementary regions of the available parameter space. 

Results are presented in this Letter in the context of the two benchmark scenarios discussed 
earlier, one in the context of NMSSM and another one in the framework of dark SUSY sce¬ 
narios. However, the search has been designed to be independent of the details of these two 
specific models, and the results can be interpreted in the context of other models predicting the 
production of the same final states. Compared to the previous version 125) , the present analy¬ 
sis has been redesigned to be sensitive to signatures with the intermediate bosons traversing a 
nonnegligible distance before decaying into a pair of muons. Such signatures can be realized 
in dark SUSY models if the mixing of the dark photon with its SM counterpart is sufficiently 
weak. In addition, the present analysis uses a dataset four times larger than the previous anal¬ 
ysis, and at a higher centre-of-mass energy, further extending the reach for signatures with 
prompt muons. 


2 The CMS detector 

This search is based on a data sample corresponding to an integrated luminosity of 20.7 fb 1 
of proton-proton collisions at a center-of-mass energy i/s = 8 TeV, recorded by the CMS de¬ 
tector in 2012. The central feature of the CMS apparatus is a superconducting solenoid of 6 m 
internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon 
pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and 
scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. Muons are measured in the pseudorapidity range \rj\ < 2.4, with detec¬ 
tion planes made using three technologies: drift tubes, cathode strip chambers, and resistive 
plate chambers. Matching muon candidates to tracks measured in the silicon tracker results in 
an accurate measurement of the transverse momentum (px). As an example, for muons with 
px < lOGeV the relative px resolution is found to be 0.8%-3.0% (depending on \rj\) and for 
muons with 20 < px < 100 GeV it is 1.3-2.0% in the barrel and better than 6% in the end- 
caps (571 . A more detailed description of the CMS detector, together with definitions of the 
coordinate system used and the relevant kinematic variables, can be found in (58l . 


3 Data selection 

The data were collected with an online trigger selecting events containing at least two muon 
candidates, one with px > 17 GeV and another with px > 8 GeV. In this analysis offline muon 
candidates are defined as particle-flow (PF) muons (57) . The PF reconstruction algorithm com¬ 
bines information from all CMS subdetectors to identify and reconstruct individual particles. 
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such as electrons, photons, hadrons or muons. 

Events are further selected by requiring at least four offline muon candidates with pj > 8 GeV 
and |?/| < 2.4 that form two oppositely charged pairs. At least one of these muons must addi¬ 
tionally satisfy the requirement of pj > 17GeV and \rj\ < 0.9, which ensures that the trigger 
efficiency is high and independent of the event topology, including effects related to overlaps 
of nearby muon trajectories. Tracks associated with a pair of opposite-charge muon candidates 
are fit for a common vertex using a Kalman filter algorithm 1591I . If the vertex is reconstructed, 
a muon pair is combined into a dimuon system if its invariant mass measured at the common 
vertex < 5 GeV and the vertex fit probability P v (p ] p') > 1%. Muon pairs failing these 

requirements are still retained for the analysis if at the point of closest approach of the two tra¬ 
jectories they are within AR(p + ,p~) = \p fl + — P},-) 2 + (<^,+ — < 0.01, where (p )L ± are 

the azimuthal angles in radians. This recovery step is designed to compensate for the reduced 
efficiency of the vertex selection for dimuons in which the two muon tracks are nearly parallel 
to each other, therefore a good efficiency is maintained for dimuon masses down to the 2m ?/ 
threshold (0.2114 GeV). For dimuons in which this is the case the point of closest approach 
is selected as the vertex position, with the additional selection requirement that the distance 
between the tracks be <0.5 mm. The dimuon kinematic variables are measured at the dimuon 
vertex position. There is no restriction on the number of ungrouped additional muons. Both 
dimuons are required to have at least one hit in the first layer of the barrel or endcaps of the 
pixel detector, and this defines an effective "fiducial" region. This requirement, along with the 
muon pj and \rj \ criteria, ensures high trigger (>96%), reconstruction, and selection efficiencies, 
with a greatly reduced dependence on the pj, //, or opening angle between the muons. 

The projected z coordinate of the dimuon system at the point of the closest approach to the 
beam line (z /(?( ) is reconstructed using the dimuon momentum. The requirement |zi ;(f , — Z 2 /(/ , | < 
1 mm is imposed to ensure that both dimuons are consistent with the same pp interaction; no 
explicit requirements are made on the impact parameter or the z coordinate at the point of 
closest approach to the beam line of the individual reconstructed muons to preserve sensitivity 
to signatures with displaced muons. 

To suppress background events in which the muons are produced in the decay of heavy quarks 
(and thus appear in jets), the dimuons are required to be isolated from other event activity 
using the criterion / sum < 2 GeV. The isolation parameter I sum is defined as the scalar sum 
of the pj of charged tracks with pj > 0.5 GeV within a cone of size A R = 0.4 centered on 
the momentum vector of the dimuon system, excluding the tracks corresponding to the two 
muon candidates. The tracks used in the calculation of / sum must also have a z coordinate at 
the point of closest approach to the beam line that lies within 1 mm of z w ,. The I sum selection 
suppresses the contamination from bb production by about a factor of 40, as estimated using a 
bb enriched control sample with one dimuon recoiling off a jet containing an unpaired muon, 
while rejecting less than 20% of events with the signal topology. 

The invariant mass m always refers to the dimuon containing a muon with pj > 17 GeV and 
|;/| < 0.9. For events with both dimuon systems containing such a muon, the assignment of 
mi w and tii 2 iiji is random for compatibility with the background modeling schema described in 
Section [5] The invariant masses of both reconstructed dimuons are required to be compatible 
within the detector resolution, specifically \nii }l}l — 7«2/jf<| < 0.13GeV + 0.065 {ni\^ + miufi)/2, 
which defines a diagonal signal region in the plane of the invariant masses of the two dimuons. 
The numerical parameters in the requirement correspond to at least five times the size of the 
core resolution in the dimuon mass. 




4 


4 Signal modeling 


Table 1: Event selection efficiencies e s im(wh 1 ,?n ai ) and t' s j iri 0?z 7i:)/ CT 7l J, as obtained from simu¬ 
lation, the geometric and kinematic acceptances a gen (wh 1 , m ai ) and a ge n(»* 7 D /CT 7D ), calculated 
using only generator-level information, and their ratios (with statistical uncertainties), for a 
few representative NMSSM and dark SUSY benchmark samples. The experimental data-to- 
simulation scale factor (Cdata / e S im, described later) is not applied. 


m hl [GeV] 
m ai [GeV] 

90 

2 

125 

0.5 

125 

3.55 

e S im [%] 

11.0 ±0.1 

21.1 ±0.1 

17.3 ±0.1 

^gen [%] 

15.9 ±0.1 

32.0 ±0.1 

26.3 ± 0.1 

^sim/^gen 

0.69 ± 0.01 

0.66 ±0.01 

0.66 ± 0.01 


m 7D [GeV] 
ct 7d [mm] 

0 

0.25 

0.5 

2 

0 

1.0 

0.5 

2 

Csim [7°] 

8.85 ±0.12 

1.76 ±0.05 

0.23 ± 0.03 

6.13 ± 0.23 

4.73 ± 0.07 

1.15 ±0.04 

^gen [%] 

14.32 ± 0.14 

2.7 ±0.06 

0.31 ± 0.03 

8.89 ± 0.28 

6.98 ± 0.09 

1.68 ±0.05 

e sim/ ^gen 

0.62 ± 0.01 

0.65 ± 0.02 

0.74 ±0.13 

0.69 ± 0.03 

0.68 ± 0.01 

0.68 ± 0.03 


4 Signal modeling 

The results from this analysis are designed to be model independent, but are also presented in 
the context of the two benchmark models introduced earlier. NMSSM simulation samples for 
benchmark model 1 are generated with PYTHIA 6.4.26 H60l . using MSSM Higgs boson produc¬ 
tion via gluon fusion gg —> H°mssm/ with the Higgs bosons decaying via H°mssm ~± 2A ( r |,, SSM . 
The masses of the MSSM bosons H°mssm an d A ( ^ 1SSM are set to the desired values for the h | 
mass and ai mass of the NMSSM bosons, respectively. The mass of H°mssm is in the range 
90 — 150 GeV (mass below 90 GeV is excluded by LEP ||37| ) and the mass of Aj^ gSM is in range 
0.25 — 3.55 GeV. Both A^ SSM bosons are forced to decay promptly to a pair of muons. Dark 
SUSY simulation samples for benchmark model 2 are generated with MadGraph 4.5.2 f6l1 
using SM Higgs boson production via gluon fusion gg -A hgM, with m^ SM = 125 GeV. The 
BRIDGE program If62l is used to force the Higgs bosons to undergo a non-SM decay to a pair 
of neutralinos, each of which decays via ni —> np ± 7 d , where m ni = 10 GeV, m„ D = 1 GeV, 
which is representative of the type of models considered l42l . Dark photons are generated with 
m 7o in the range 0.25-2.0 GeV and a decay length ct 7d in the range of 0-20 mm. Each of the two 
dark photons are forced to decay to two muons, while both dark neutralinos escape detection. 
The narrow-width approximation is imposed by setting the widths of the dark photons to a 
small value (1CU 3 GeV). 

All benchmark samples are generated using the leading-order CTEQ6.6 Ii63l set of parton dis¬ 
tribution functions (PDF), and are interfaced with PYTHIA using the Z2* tune 11641 for the un¬ 
derlying event activity at the LHC and to simulate jet fragmentation. 

The signal samples are processed through a detailed simulation of the CMS detector based on 
Geant 4 l65l and are reconstructed with the same algorithms used for data. Table [l] shows 
the event selection efficiencies t' s j m obtained using the simulated benchmark samples for a few 
representative choices of (m^, m a ,), and (m 7D , ct 7d ). To provide a simple recipe for future rein¬ 
terpretations of the results in the context of other models, the variable a gen is separately defined 
as the geometric and kinematic acceptance of this analysis calculated using only generator-level 
information. It is defined by selecting events containing at least four muons with pj > 8 GeV 
and \ f] | < 2.4, with at least one of these muons having pj > 17 GeV and \ f] < 0.9. The new 
light boson must also decay with transverse decay length L xy < 4.4 cm and longitudinal decay 
length L z < 34.5 cm (both defined in the detector reference frame), to satisfy the "fiducial" 












5 


region of the analysis. Table [I] shows a gen along with the ratio e' s j m /ct gen . 


5 Background estimation 

The SM background for this search is dominated by bb production and has small contributions 
from the electroweak production of four muons and direct J /xp pair production. The leading 
part of the bb contribution is due to b quark decays that result in a pair of muons, via either 
the semileptonic decays of both the b quark and the resulting c quark, or via resonances, i.e. 
co, p, (p, J/xp. A smaller contribution comes from events with one genuine dimuon candidate 
and a second dimuon candidate containing one muon from a semileptonic b quark decay and 
a charged hadron misidentified as another muon. 

Using data control samples, the bb background is modeled as a two-dimensional (2D) tem¬ 
plate in the plane of the invariant masses of the two dimuons. The tem¬ 

plate describing the 2D probability density function is constructed as a Cartesian product 
Bi 7 (nii m ) Bs(m 2 fiji), where the £>17 and Bg templates model the invariant mass distributions 
for dimuons with and without the requirement that the dimuon contains at least one muon sat¬ 
isfying pj > 17GeV and \tj\ <0.9 respectively. The Bn shape is measured using a data sample 
enriched with bb events containing exactly one dimuon and one additional muon, under the 
assumption that the decay of one of the b quarks results in a dimuon pair containing at least 
one muon with pj > 17 GeV and \i]\ < 0.9, while the other b quark decays semileptonically 
resulting in the additional muon with pj > 8 GeV. For the By shape, a similar sample and 
procedure is used but the dimuon is required to have both prongs with pj > 8 GeV, while 
the additional muon must have pj > 17GeV and \rj\ < 0.9. The two templates are required 
as the shape of the dimuon invariant mass distribution depends on the pj thresholds used to 
select the muons and whether the muons are restricted to the central (\tj\ < 0.9) region or 
can be in the full acceptance range (|?/| < 2.4), as a result of the differences in the momentum 
resolution of the barrel and endcap regions of the tracker. The Bn and By distributions are 
fitted with a parametric analytical function using a sum of Bernstein polynomials and Crystal 
Ball functions 1661 describing resonances. These event samples do not overlap with the sample 
containing two dimuons that is used for the main analysis and they have negligible contribu¬ 
tions from non-bb backgrounds. Once the B bb {m^ }l}l , m^) template is constructed, it is used 
to provide a description of the bb background shape in the signal region. This technique as¬ 
sumes that each b quark fragments independently and that if the shapes of these distributions 
are measured using data samples with kinematics very similar to that of the background events 
then the effects of residual kinematically-induced correlations are small (albeit weakly induced, 
the shape depends on the b jet pj and the pj of the two b jets in background bb events tend 
to be similar). The background template is validated in a region where both dimuons fail the 
hum < 2 GeV requirement and good agreement with data is observed. 

The data events that satisfy all analysis selections but fail the ni \ w ~ 1112 , 1 ,, requirement are used 
to normalize the B^Jini^, m 2 }iyi) template. This selection yields nine events in the off-diagonal 
sideband region of the (ffli w ,m 2 „ f( ) plane, leading in the diagonal signal region to an expected 
rate of bb background events of 2.0 ± 0.7. This is essentially (9 ± y/9) x 0.18/0.82, where 0.18 
and 0.82 correspond to the integral of the areas under the background template inside and out¬ 
side the signal diagonal region, respectively. These nine events in the off-diagonal sidebands 
of the (ffli w , iiiiiiu) plane are shown as white circles in Fig. [I] 

The contribution from direct J/xp pair production is estimated using another data control sam¬ 
ple. Events are selected with a trigger that requires at least three muon candidates, two of which 
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5 Background estimation 
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Figure 1: Distribution of the invariant masses mi (( „ vs. m 2 ^ for the isolated dimuon events 
following the application of all constraints except the ~ requirement of compatibility 

within the detector resolution. The compatible diagonal signal region (outlined with dashed 
lines) contains one data event (triangle) at mi MM = 0.33 GeV and = 0.22 GeV. There are 
also nine data events (white circles) which fail the m-\^ — compatibility requirement. The 
color scale indicates the expected SM background in range 2 m fl < m\ m , ni 2 ^ < 2 m T . 

have a common vertex and an invariant mass consistent with that of the J/ip particle. Events 
are further required to contain at least four reconstructed muons with pj > 3.5 GeV, which 
form dimuon pairs. This control sample does not specifically require that the dimuons satisfy 
the requirement f sum < 2 GeV since I SU m is used to separate the contribution of "prompt" and 
"nonprompt" (from b quark decays) J/ip in data. Finally, both dimuons are required to have an 
invariant mass between 2.8 and 3.3 GeV. Following these requirements the data sample con¬ 
sists of events containing prompt and nonprompt J/ip. To subtract the nonprompt component, 
two independent methods have been studied: the first one divides the control sample based on 
the values of the isolation variable / sum for each of the two dimuons in each event. The number 
of events in which both dimuons satisfy the requirement / sum < 2 GeV is extrapolated from 
the regions in which at least one of the dimuons fails this requirement. The second approach 
uses the lifetime of ]/ip candidate, calculated under the hypothesis of it being produced at the 
beam line, as a discriminating variable. The data distribution is fitted in the isolated region 
using prompt and nonprompt templates from simulation and nonisolated sideband in data, re¬ 
spectively. Both approaches give consistent results within the associated uncertainties and the 
results of the isolation-based method are used in the final analysis. There are two mechanisms 
for the production of prompt double j/ip events: single- and double-parton scattering (SPS and 
DPS, respectively), corresponding to whether the two J /ip mesons are produced from one or 
two independent parton interactions. The number of prompt events in the control region are 
further separated into SPS and DPS components using the J /ip rapidity difference as the dis- 
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criminating variable. Finally, the data-to-simulation normalization factor and the fraction of 
SPS and DPS events are extrapolated from the control to the signal region, resulting in a final 
estimation for the contribution from prompt double J/ip events of 0.06 ± 0.03 events. 

The contribution from other SM processes (low mass Drell-Yan production and pp —» Z/ 7 * —t 
4/<) is estimated with the CalcHEP 3.6.18 generator Ii67j using the HEPMDB infrastructure If 68 l . 
and is found to be 0.15 ± 0.03 events in the entire signal region. The combined expected back¬ 
ground contribution to the diagonal signal region is 2.2 ± 0.7 events. This background is repre¬ 
sented by the color scale in Fig. |T] 

6 Systematic uncertainties 

The selection efficiencies for the offline muon reconstruction, trigger, and dimuon isolation 
requirements are obtained from simulation, and are corrected with scale factors derived from 
comparison between data and simulation using Z —> //// and J/ip —» ftp samples. The scale 
factor per event is found to be e^ata /%im = 0.93 ± 0.07 and it accounts for the differences in the 
efficiency of the trigger, the efficiency of the muon reconstruction and identification for each 
of the four muon candidates, and the combined efficiency of the isolation requirement for the 
two dimuon candidates. The estimate accounts for correlations associated with the presence of 
multiple muons per event. The main systematic uncertainty is the offline muon reconstruction 
(4.1%), which includes an uncertainty (1% per muon) to cover variations of the scale factor as 
a function of the muon pj and if. Other systematic uncertainties include: the uncertainty in 
dimuon reconstruction effects related to overlaps of muon trajectories in the tracker and in the 
muon system (3.5%), the trigger efficiency (1.5%), the uncertainty in the efficiency caused by the 
modeling of the tails in the dimuon invariant mass distribution that arises from the requirement 
that the two dimuon masses are compatible (1.5%), and the dimuon isolation (negligible). The 
uncertainty in the integrated luminosity of the data sample (2.6%) l69l is also included. All 
uncertainties quoted above are related to variations in the signal efficiency due to experimental 
selection and sum up to 6.3%. The uncertainties related to variations in the signal acceptance 
due to the model include: the uncertainties related to the PDFs and the knowledge of the strong 
coupling constant oc s , which are estimated by comparing the PDFs in CTEQ 6.6 ||63l with those 
in NNPDF2.0 Il70l and MSTW2008 [|7TH , following the PDF4EHC recommendations If72l [73} . 
Using the analysis benchmark samples, they are found to be 3% for the signal acceptance. The 
variation of the renormalization and factorization scales has a negligible effect. In addition 
a re-weighting procedure is applied to the Higgs boson pj spectra in the benchmark signal 
samples to reproduce the NNFO+NNEF prediction ||74| and account for possible changes to 
the analysis acceptance. Only a weak sensitivity to this is expected and the result of the re¬ 
weighting procedure is limited by the statistical uncertainty of the simulated samples ( 2 %), 
which is therefore included as a conservative systematic uncertainty on this effect. Thus, the 
total systematic uncertainty in the signal acceptance and selection efficiency is 7.3%. 


7 Results 

After the full analysis selection is applied to the data sample, one event is observed in the 
diagonal signal region, as shown in Fig. [lj This is consistent with the expected background 
contribution of 2.2T0.7 events. 

For future reinterpretations of this analysis, the results can be presented as a 95% confidence 
level (CF) upper limit on a (pp —> 2a + X) £> 2 (a —> 2p) a gen = N^)/ (£f), where a gen is the 
generator-level kinematic and geometric acceptance defined earlier. The calculation uses the in- 
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7 Results 


tegrated luminosity C = 20.7fb 1 and central value f of the ratio r = e^ata/ a gen = 0.63 ± 0.07. 
The ratio includes a scale factor correcting for experimental effects not included in the simula¬ 
tion and its uncertainty covers the variation in the ratio over all the benchmark model points 
used. The limit is calculated as a function of the dimuon mass using the CLs approach 117511761. 
The chosen test statistic is based on the profile likelihood ratio and is used to determine how 
signal- or background-like the data are. Systematic uncertainties are incorporated in the anal¬ 
ysis via nuisance parameters with a log-normal probability density function and are treated 
according to the frequentist paradigm. The overall statistical methodology used in this anal¬ 
ysis was developed by the ATLAS and CMS Collaborations in the context of the LHC Higgs 
Combination Group and is described in |3]. [7Zl- The obtained limit as a function of dimuon 
mass m^m can be conveniently approximated as a constant everywhere except the vicinity of 
the observed event, where it follows a Gaussian distribution: 


N(m }l}l ) < 3.1 + 1.2 exp 


( {m m — 0.32) 2 \ 
^ 2 x 0.03 2 )' 


resulting in 


cr (pp —> 2a + X) B 2 {a —¥ 2 p) a gen < 0.24 + 0.09 exp 


( ( m uu — 0-32) 2 \ 
^ 2 x 0.03 2 J ’ 


where is measured in GeV and the cross-section limit is expressed in femtobarns. This 
limit is applicable to models with two pairs of muons coming from light bosons of the same 
type with a mass in the range 2m ^ < rn a < 2 m T , where the new light bosons are typically 
isolated and spatially separated (so as to satisfy the isolation requirements). 

The weak model dependence of the ratio r allows for a simple reinterpretation of the results in 
other models. This requires calculating a gen , as defined earlier, and then the full efficiency e^ata 
can be calculated by multiplying a gen by the ratio r. 

There are certain subtleties that must be taken into account when reinterpreting the model- 
independent results of this analysis in the context of other models, particularly with the isola¬ 
tion requirement. An event should be considered to satisfy the selection requirements if there 
are at least two well isolated 7 d decaying to muon pairs. Experimentally, isolation is based on 
charged tracks but it may be insufficient to just require the absence of generator-level charged 
particles in the isolation cone. For example a neutral pion decaying to a pair of photons, that 
convert into electrons, may result in the reconstruction of one or more tracks. This would be 
particularly relevant for models with more than two dark photons produced in the same event, 
some of which may decay to hadrons or electrons. In this case the safest approach is to require 
that there are no particles with pj > 0.5 GeV within the 7 d isolation cone. This restriction 
would result in a more conservative limit but it would be robust against these effects. 

The results from this analysis are also interpreted in the context of the NMSSM and the dark 
SUSY benchmark models, and 95% CL upper limits on the product of the cross section and 
branching fraction are derived. In these models both the Higgs boson production cross section 
and the branching fractions can vary significantly, depending on the choice of parameters. In 
the absence of broadly accepted benchmark scenarios, the production cross sections in these 
examples are normalized to that of the SM Higgs boson with a mass of 125 GeV I78l . 

In the case of the NMSSM benchmark scenario, the production cross sections and branching 
fractions for hi and h 2 can vary substantially depending on the chosen parameters. An ex¬ 
act interpretation of these results requires evaluating the experimental acceptance using the 
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Figure 2: Left for benchmark model 1: 95% CL upper limits from this search for the NMSSM 
scenarios with m 3] = 0.25 GeV (dashed curve), m a , = 2 GeV (dash-dotted curve) and m ai = 
3.55 GeV (dotted curve) on u(pp —> \\\/ 2 —> 2a\) B 2 {a\ —> 2}i) as a function of in the range 
86 < ?77hj < 125 GeV and of «%, for m h , > 125 GeV. As an illustration, the limits are compared 
to the predicted rate (solid curve) obtained using a simplified scenario with cr(pp —> h, -4 
2a i) — 0.008 UgM, which yields predictions for the rates of dimuon pair events comparable to 
the obtained experimental limits, and B{a\ —> 2 }i) = 7.7%. The chosen B(a\ —> 2}i) is taken 
from IfTTH for m ai = 2 GeV and tan /3 = 20. Right for benchmark model 2: 95% CL upper 
limits (black solid curves) from this search on u(pp —> h —> 27 D + X) B(h —> 27 ^ + X) (with 
7 /z ni = 10 GeV, 7 /z no = lGeV) in the plane of two of the parameters (e and ;h 7d ) for the dark 
SUSY scenarios, along with constraints from other experiments l42l - |56l showing the 90% CL 
exclusion contours. The colored contours represent different values of £>(h —>■ 270 + X) in the 
range 0.1-40%. 


generator-level acceptance for each of the two hi 2 bosons, and then using the measured upper 
limit on the sum of two contributions to derive limits for any choice of NMSSM model param¬ 
eters. To present results in a fashion allowing for straightforward interpretation, we note that if 
one of the two CP-even Higgs bosons is the 125 GeV state observed at the LHC, then the other 
one is either lighter or heavier. In the NMSSM it is typical that one of the two has approxi¬ 
mately the SM production cross section and a small £?(h, —> 2a), whereas the other one has 
a suppressed production rate and large B(hj —r 2a) due to its large singlet fraction. In Fig. [ 2 ] 
(left) the limit at each mass point is calculated taking the CP-even Higgs boson with the corre¬ 
sponding mass as the only source of signal events; the curve below 125 GeV applies to NMSSM 
models in which < m^ 2 = 125 GeV, with hi decays dominating the rate of 4 j.i events. The 
limit at 77 Z;, = 125 GeV corresponds to the case where 125 GeV = < m h 2 , with hi decays still 

responsible for the vast majority of signal-like events. The points above 125 GeV correspond to 
model points for which only I 12 (77Zh 2 > z«h = 125 GeV) is allowed to have a sizeable rate of 
observable 4 }i events. Finally, for models with m^ 2 > 150 GeV, the limit at 150 GeV can be used 
as a conservative estimate of the production rate limit. In each of these scenarios it is possible 
that the other Higgs boson also contributes some fraction of the 4w signal events, in which case 
the limit shown is more conservative than would be given by an exact evaluation. 

In the case of the dark SUSY scenario, a 95% CL limit on the product of the Higgs boson 
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8 Summary 


production cross section and the branching fractions of the Higgs boson (cascade) decay to 
a pair of dark photons is determined. The limit set in the (m 7D , e) plane from this analysis is 
shown in Fig. [2] (right), along with limits from other experimental searches, where the lifetime 
is directly related to the kinetic mixing parameter e and the mass of the dark photon m lD via 
r 7 D (e, m 7o ) = e~ 2 where /(m 7l J is a function that depends only on the mass of the 

dark photon H79| . The significant vertical structures in the limits visible in Fig. [2] (right) arise 
because the total width of the dark photon varies rapidly in those mass regions due to reso¬ 
nant decays to hadrons. This search constrains a large, previously unconstrained area of the 
parameter space. Unlike the other results in the figure, the CMS and ATLAS limits are model- 
dependent and only valid under the assumption that 13(h —> 2n\ —» 4// + X) 7 ^ 0. The recent 
ATLAS analysis 11421 focused on highly displaced objects and these searches therefore probe 
different regions of the available parameter space. 


8 Summary 

A search for pairs of new light bosons produced in the decay of a Higgs boson, that subse¬ 
quently decay to pairs of oppositely charged muons, is presented. One event is observed in the 
signal region, with 2.2 ± 0.7 events expected from the SM backgrounds. A model independent 
upper limit at 95% CL on the product of the cross section, branching fraction, and acceptance 
is obtained. This limit is valid for light boson masses in the range 2m,, < m a < 2m r . The 
obtained results allow a straightforward interpretation within a broad range of physics models 
that predict the same type of signature. The results are compared with two benchmark models 
in the context of the NMSSM and dark SUSY, including scenarios predicting a nonnegligible 
light boson lifetime. 
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